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Abstract 



The probability distribution of lens image separations is calculated for 
the "standard" gravitational lensing statistics model in an arbitrary, flat 
Robertson- Walker universe, where lensing galaxies are singular isothermal 
spheres that follow the Schechter luminosity function. In a flat universe, the 
probability distribution is independent of the source distribution in space and 
in brightness. The distribution is compared with observed multiple-image lens 
cases through Monte-Carlo simulations and the Kolmogorov-Smirnov test. 
The predicted distribution depends on the shape of the angular selection bias 
used, which varies for different observations. The test result also depends 
on which lens systems are included in the samples. We mainly include the 
lens systems where a single galaxy is responsible for lensing. We find that 
the "standard" model predicts a distribution which is different from the ob- 
served one. However, the statistical significance of the discrepancy is not large 
enough to invalidate the "standard" model with high confidence. Only the 
radio data reject the model at the 95% confidence level, which is based on 
four/three lens cases. Therefore, we cannot say that the observational data 
reject the "standard" model with enough statistical confidence. However, if 
we take the velocity dispersion of dark matter without the conversion factor 
(3/2) 1 / 2 from that of luminous matter, the discrepancy is quite severe, and 
even the ground-based optical survey data reject the "standard" model with 
^ 90% confidence. 

PACS numbers: 95.30.S, 98.62.S, 98.54.Aj, 98.65.Dx, 98.80.Es 
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I. INTRODUCTION 



Determining the geometry of the universe and the distributions of masses therein has 
always been one of the most important goals of astronomy. Even before the actual discov- 
ery of gravitational lenses, it was known that lens systems could provide us much useful 
information on these fundamental questions 0. 

As more high-redshift multiply imaged QSOs were discovered, Turner et al. [|l| calculated 
the statistical properties of multiple-image gravitational lenses in the standard Robertson- 
Walker cosmology, especially the probability of a QSO being multiply imaged and the mean 
angular separation of images [IJ. They compared the predicted statistical properties of the 
lens systems with observations and concluded that although more small-separation 1") 
lenses are expected, the distribution of the image separation, with the effect of resolution 
bias considered, is compatible with the observed lens systems. Gott et al. extended this 
work to an arbitrary Robert son- Walker universe || . The probability of lensing and the mean 
image separation were expressed as functions of the source redshift. However, the explicit 
value of the mean image separation was not calculated, and only the functional form was 
compared with the data to get the limits on the cosmological parameters. One interesting 
fact recognized in this work, and partly in Ref. H, is that in a flat universe filled with 
singular isothermal sphere (SIS) galaxies, the mean separation of images is independent of 
the source distance (or redshift), regardless of the value of the dimensionless cosmological 
constant, Ao = A/3Hq where H is the Hubble constant. 

Turner and Fukugita et al. |J subsequently realized that the probability of lensing, 
rather than the image separation, is more sensitive to Ao- In a series of papers they 
estimated how many multiple-image lens systems should be discovered for given QSO sam- 
ples in a Ao-dominated universe and found the present number of observed lens systems 
exclude Ao ^ 0.9. This limit is considered as one of the few observational limits on the 
elusive cosmological constant. 

However, if we naively use the various galaxy parameters of Ref. || to recalculate the ex- 
pected mean separation in a flat universe following Ref. [Q, we get 1.5". This is significantly 
smaller than the observed image separations in known single-galaxy lens systems (systems 
in which a single galaxy is responsible for lensing), e.g., PG1115, Q2237, Q0142, H1413, 
MG0414, and B1422 (Table |) whose mean value is 1.95"and whose dispersion 0.70". If we 
take the most favorable galaxy parameters allowed within the uncertainty interval in Ref. 
@, the mean separation goes up to 1.7", still somewhat small compared to the observed 
image separations. However, the angular resolution bias, larger splitting lens systems being 
easier to find, which was not considered in Ref. Q and which was minimally incorporated 
in Refs. J7j| and [§] might change this interesting discrepancy. Also, a better statistical test 
is needed because of the non-standard shape of the distribution and the small sample size. 

The usual criticism against utilizing image separation data has been its dependence on 
a combination of cosmological parameters |4j]. However, this can be an advantage if we 
want to test the curvature of the universe or the consistency of the lensing statistics model 
without being affected by the uncertainties in the cosmological parameters. Moreover, the 
distribution of the image separations is relatively independent of the magnification bias 
which introduces major uncertainty when the probability of lensing is used as a cosmological 
diagnostic tool. 
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In this work, we adopt the same assumptions and parameters for galaxies as in Ref. 
@ to get the probability distribution of the image separations, incorporating all possible 
angular resolution biases present in various lens surveys. By comparing this distribution 
with observed single-galaxy lens data (Table 1), we test the "standard" gravitational lensing 
statistics model against the observations. We choose a flat universe with arbitrary Ao as the 
background universe. 



II. PROBABILITY DISTRIBUTION OF IMAGE SEPARATION 

A flat universe with a non-zero cosmological constant can be described by the Robertson- 
Walker metric 01 



ds 2 



fdP + ^rWx 2 + a 2 oX 2 (d6 2 + sin 2 6d<j> 2 )]. 
a 



(1) 



The expansion factor a(t) has as its present value a = Rq/A, where Rq = c/H and 
A = (\Q + Ao — l)) 1 ^ 2 if the value inside the square-root is not zero and A = 1 if it is zero. 
Among several relevant cosmological distances, we use the parametric distance x, which is 
equivalent to the (filled-beam) standard angular distance (see Ref. ||). 

Galaxies are modeled as SIS; therefore, their lensing properties are determined by a 
one-dimensional velocity dispersion, a, which is related to the galaxy luminosity by L/L* = 
(cr/cr*) 7 where 7 = 4 for E and SO galaxies and 7 = 2 for for S galaxies The luminosities 
of the galaxies follow the Schechter function 



<&(L)dL = <F 



L* 



exp(-L/L* 



dL 

'I* 



(2) 



where a g = —1.1 |U|]. The exact values of the typical luminosity of a galaxy, L*, and its 
normalization constant, $*, are important in lensing probability estimates, but are irrelevant 
in this work because only the functional shape of the distribution matters. If there are galaxy 
evolutions, $* and L* become functions of x ( or redshift). 

Under these assumptions, the lensing probability by one type of galaxy is simply given 
as § 



dV = 7rag(a*) 2 $* - 



-L/L* 



(Xs ~ Xl ) 2 X 2 
Xs 



d \j7) d XL 



(3) 



where \s an d Xl are the parametric distance to the source and the lens, respectively, and 
a* = A7r(a*/c) 2 , the bending angle for an L* galaxy. Since the image separation A# = 
2a*(L/L*) 2 ^(xs — Xl)/xs, the probability can be rewritten as a function of AO or ip = 
A6/(2a*): 



d 2 r = | 7 a3(a*) 2 $V^ +1 ^ 2+1 exp 



( XS-XL j- 7 /2 7/2 

Xs 



X 



(XS ~ Xl) 2 x\ 



Xs 



Xs ~ XL 

Xs 



- 7 K+i+4/ 7 )/2 



dipdxL- 



(4) 
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The total probability of lensing by all types of galaxies is 



T 



E 

=E,S0,S 



xs d 2 r t 
d<Pi / d XL-i—; — 
o Jo dipidXL 



(5) 



where G?Tj and the other symbols with subscript i refer to the respective values for E, SO, 
and S galaxies. 

All lens observations and surveys suffer from some kind of angular selection bias mainly 
due to the finite angular resolution. If this bias is taken account, the probability distribution 
function (PDF) of observed image separations becomes 



dP 
dA6 



1 rxs 

T JO 



E 



d 2 - 



dxi 



dxLdA9 / 

E^q?)- V?[r(a g + 1 + 4 7 r 1 )]~ 1 f xs dxLGjjxL, fi) 
£i 7^ir(« 9 + 1 + 47," 1 )]- 1 /o°° % Io Xs dXLGiixL, <*) 



(6) 
(7) 



with 



Gi(XL,<Pi 



f s (A9) 
x exp 



(Xs ~ Xl)Xl 



-( 



Xs 
Xs 



Xs 



Xs ~ XL 



7i (a s +l)/2+l 



Xs ~ XL 



7i/2 7i /2 



(8) 



where T is the gamma function, f s (A8) is a function describing the angular resolution bias 
with ipi = A6>/(2a*), and F* = ir§*(a*) 2 Rlr(a g + 1 + ^j' 1 ) is the dimensionless parameter 
for the effectiveness of matter in producing double images for each type of galaxy |T|[7|. 
We take the values for the various parameters from Ref. ||. The cumulative distribution 
function (CDF) 



P(A6) 




dA6' 



(9) 



is calculated from dP/dAO. 

It can be proved with ease that in a flat universe the distribution is independent of xs 
(or zs, the redshift of the source), regardless of the value of Ao or the functional form of the 
angular resolution bias {§. This implies, regardless of the source distribution in space, any 
multiple-image lensing event should have the same image separation distribution as long as 
the universe is flat, which makes various statistical tests very simple. However, this is true 
only when there is no evolution in the lensing galaxies. 



III. OBSERVED MULTIPLE-IMAGE SYSTEMS AND ANGULAR SELECTION 

BIAS 

The image separation is not affected by the absolute brightness of the sources as is the 
case with the lensing probability. However, it is affected by the finite angular resolution 
and by the dynamic range of the real observations; therefore, the angular selection bias 
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is a function of the image separation and the relative brightness of each image. Hence, 
to analyze statistically the distribution of image separations of lens systems discovered in 
specific surveys, one should carefully consider both effects as in Ref. | TT| . However, since 



we approximate the lens by SIS in this work and since a given SIS lens produces always the 
same image separation regardless of the source position, the angular selection bias will be a 
function only of the image separation. 

Since most of the lens cases are not discovered in one systematic survey and since each 
survey or observation is done under different conditions, a single universal function does not 
exist which can faithfully represent the angular selection effect in various different observa- 
tions. Hence, we must use different specific angular selection biases for different observa- 
tions. In this work, we consider four kinds of selection biases suitable for specific surveys: 
no selection bias (NO), the one for the HST snapshot survey |Tl|JT^] (HST), four similar 
selection biases appropriate for ground optical lens surveys [0 (PSF1.0, PSF0.7, EYE1.0, 
EYE0.7) one for a radio survey [15,3] (MG-C). Although the original forms for the selec- 
tion biases HST, PSF1.0, PSF0.7, EYE1.0, and EYE0.7 depend on the brightness ratio of 
the images, we can disregard this dependence for the reasons explained above. Hence, we 
use Kochanek's completeness function [|12j as the selection bias function because the image 
separation statistics are independent of the source brightness or the redshift distributions. 

The lens samples and their corresponding selection biases with descriptions are listed in 
Table |H| and shown in Fig. [I]: Inverted big triangles, small triangles, and circles denote the 
image separations of systems with a single galaxy lens or unknown lens, systems with a lens 
other than a single galaxy, and radio rings. We use only the accepted cases of multiple-image 



systems according to the criterion of Ref. ||16||. Also, we basically exclude systems like 0957 



(lensed by a galaxy and a cluster) and MG2016 (lensed by two galaxies) because we only 
consider single-galaxy lensing. However, the two systems are included in a special case for 
comparison. In systems like H1413 and MG0414, we do not know what kind of object is 
responsible for lensing. However, excluding those systems could introduce a bias against 
a small-separation system which has a less massive galaxy and, therefore, less chance of 
the lens being detected. For that reason, we include them as single-galaxy lens systems. 
Although system Q1208 is classified only as a proposed case in Ref. it is meaningless 



to discuss the HST snapshot survey result without it |TT|, and we include it only for the 
HST survey. The test result regarding the HST survey should be taken with some cautions. 
Radio rings are included only in MG-C radio survey samples. 



IV. STATISTICAL TEST 

Figure [I] shows the PDF's of the image separations for each angular selection bias. All of 
them peak around sub-arcsecond separations. However, all of the observed single-galaxy lens 
systems (big triangles) have image separations larger than 1", showing a discrepancy with 
theoretical expectations. To access the degree of discrepancy statistically, we use Monte- 
Carlo (MC) simulations [|TTJ and the Kolmogorov-Smirnov (KS) test. 
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A. Monte-Carlo Test 



We define a likelihood function [IT1 



In 



dP 



(10) 



and generate 10 4 sets of four to six random A#j's following the calculated PDF and their 
respective i values. These are compared with the I values for the observed A0j's, and the 
number of MC sets having i values less than that of the observed sample is counted. The 



result is expressed in (absolute) probability in Table III. We can see that most models are 
not rejected and that only the no-selection-bias model is rejected if we include the large 
separation, non-single-galaxy lens cases, MG2016 and Q0957. Also, the model with the 
MG-C selection bias is rejected at ~ 90% level by the MG-C band survey samples. 



B. Kolmogorov-Smirnov Test 

The MC simulation is quite versatile, being applicable even to the case where the dis- 
tribution differs for each sample data. However, the test mainly focuses on whether the 
data occur at the probable part of the distribution. It does not test whether the data are 
clustered or not. The KS test is complementary to the MC simulation because it has good 
discriminative power on the clustering of the data. Hence, we use the KS two-sided test 
to find out if the calculated distribution is significantly different from that of the observed 
samples. This is possible because in a flat universe the distribution of separation angles is 
the same regardless of the source redshift, thus making the KS test easily applicable. 

The CDF's constructed from the observed lens samples (jagged lines) and from the 
calculated distribution (smooth curves) are shown in Fig. |^: (a) the dot- dashed curve is 
CDF for no selection bias (NO), the solid line CDF from lens systems O-ALL, the dotted 
line from O-SINGLE, and the dashed line from O/R-SINGLE; (b) the dotted curve for bias 
HST and the solidline from corresponding lens systems; (c) the dotted curve for bias PSF1.0, 
the short- dashed curve for bias PSF0.7, the long- dashed curve for bias EYE1.0, the dot-dashed 
curve for bias EYE0.7 and the solid line from corresponding lens systems; (d) the dotted 
curve for bias MG-C and the solid line from corresponding lens systems. The comparisons 
between the observed and the calculated CDF's show the largest differences for the NO bias 
and the MG-C bias models. 

The results of the KS test applied to each bias case are shown in Table [IV|. The numbers 
show the confidence of rejection, and only NO-O-ALL and MG-C are rejected with more 
than 95% confidence. This is in agreement with the MC test results. Hence, we can say 
that although the observed optical lens systems show somewhat different distribution from 
theoretical ones, the difference is not statistically significant. Only the radio lens systems 
show statistically significant inconsistencies with the "standard" lens statistics calculation. 
However, this is based on only four systems. If we dismiss lens system MG2016 which has two 
lensing galaxies, the sample consists of only three systems. Hence, although a statistically 
significant conclusion can be drawn formally, it is based on a few lens systems and should 
be taken with cautions. 
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C. Without the (3/2) 1 / 2 Factor 



Turner et al. JTJ] were the first to introduce the correction factor (3/2) x l 2 to convert 
the observed velocity dispersion (of luminous matter) into that of dark matter for E and 
SO galaxies (also see Ref. Q). However, Kochanek argues against this correction and 
estimates a "90% confidence range" in the velocity dispersion which is significantly smaller 
than the (3/2) 1//2 corrected value. Hence, we repeat the calculation and the test without 
the (3/2) 1 ' 2 factor. Since the image separation in a SIS lens is directly proportional to the 
square of the velocity dispersion, the mean image separation get smaller by 67% and the 
whole distribution moves to smaller values. The PDF's with and without the (3/2) 1//2 factor 
are compared in Fig. |3] (solid for with and dotted for without). The result of the MC test 
is summarized in Table |TJ: All models, except HST and EYE1.0, have probabilities less 
than 10%, and the no-selection-bias models and MG-C model have probabilities less than 
5%. The result of the KS test is shown in Table [TV]: Again, the confidence of rejection is 
^ 90% for all models except HST, and no-selection-bias models and the MG-C model are 
rejected by the observed lens systems with more than 97% confidence. To summarize, if 
the (3/2) 1//2 factor is really unnecessary, the observed lens data reject the "standard" lens 
statistics models, marginally or strongly depending on the angular selection biases and the 
corresponding data sets used. 



V. DISCUSSION 

There are several factors which can contribute to the uncertainties in the probability 
distribution and in the statistical test. Firstly, the test depends sensitively on the observation 
sample due to the small number of appropriate lensing cases. From the analysis, we have 
seen that the conclusion does depend on which observed lens systems are chosen. This 
problem will disappear when we have enough lens systems. 

We may increase the observed sample size by including other lens systems like Q0957+561 
where the lensing galaxy is aided by a cluster. However, then the theoretical calculation 
should be modified accordingly to incorporate such multiple deflector cases. We think it 
is misleading to include lens system like Q0957 in the observed sample and to compare it 
with the theoretical distribution calculated by assuming only one galaxy as the deflector. 



However, this was done in many previous works p|Jll|JT! 

It is also possible to use the measured redshifts of the deflectors and to calculate the 
distribution of image separations for a given lens redshift and source redshift. This approach 
would be useful if we had enough lens systems because the distribution in this case would 
depend on flo and Ao individually and would be relatively insensitive to magnification bias. 
However, the number of lens systems found so far is too small to make any statistically 
significant discrimination We need more clean lens cases to make the test meaningful. 

Secondly, the presence of cores in the galaxies and asphericity certainly affects the dis- 
tribution. The core radii of galaxies are generally small enough to make SIS a good ap- 



proximation P,ni|. If the core radii are not smaller than the typical image separations in 



the lens plane, the expected separation of images will be smaller and the discrepancy will 
be larger. The ellipticity can affect lensing frequencies via the magnification bias [|H|,pf. 
However, the distribution of image separations is relatively insensitive to the magnification 
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bias — the magnification bias can affect the distribution only indirectly through the angular 



selection bias which is weakly linked to the magnification bias [ 12| . It is possible to have 



high magnification events having much smaller separations than the characteristic values of 
the SIS cases. However, these types of events are expected to be rare |7j . 

The rather strong disagreement between the observation and the theoretical expectation, 
when we omit the (3/2) x l 2 factor, is not apparent in the extensive maximum likelihood study 
of Kochanek [F2| . The main reason for this difference is the observed lens systems used. Lens 
systems Q1208+101, B1938+666, and B0218+356, all having <^ 1" separations, are included 
in Kochanek's sample — they are omitted from our sample because they do not pass the 
criterion for the accepted case or because the image separation is not known yet. 

Now, we can think of possible reasons for the discrepancy between the probability distri- 
bution from the theory and that from the observational data. The main ingredients of the 
"standard" model can be divided into two parts: one regarding cosmology and the other re- 
garding galaxies. To determine if different cosmological model would change the test result, 
we tried two non-flat cosmological models (the open model with ft = 0.1, A = and the 
closed model with Q — 2, X — 0). We applied a MC simulation test to these models. The 
probabilities from the MC simulations increased by 10 ~ 30% for the closed universe and 
decreased by 20 ~ 30% for the open universe because the image separation is smaller in an 
open universe and larger in a closed universe than in a flat universe [Q. Yet, the result was 
inconclusive again. Hence, the cause of the inconsistency between the "standard" model 
and the observation is likely to be in our poor understanding of galaxies rather than in the 
cosmological models. If this is the case, all previous works using the "standard" lensing 
statistics model to probe the cosmological parameters should be taken with caution. 



VI. SUMMARY 

We calculated the distribution of image separations in the "standard" gravitational lens 
statistics model where the universe is flat and lensing galaxies are modeled by SIS's which 
follow the Schechter luminosity function and whose comoving density is kept constant. The 
calculated distribution of image separations, incorporating the angular selection bias, is 
compared, through the MC and the KS tests, with that seen in the observed lens systems. 

The calculated distribution implies that a significant fraction of observed single-galaxy 
lens systems should have image separations <^ 1", whereas most of the observed lens sys- 
tems have separations larger than 1". However, the distribution is wide enough that the 
observations do not reject the "standard" lensing statistics model with enough statistical 
significance. The model without consideration of the angular selection bias is nearly ruled 
out, and the radio data strongly rules out the "standard" model. However, this conclusion, 
although statistically significant, is based on only four or three lens cases. 

If the correction factor (3/2) 1//2 — for converting the observed velocity dispersion to that of 
dark matter — is not used in the lens statistics calculations, the "standard" model is mostly 
rejected by the observations, although the ground-based survey is marginally compatible 
with the "standard" model. Also, larger core radii produce smaller image separations, 
and the discrepancy between the "standard" model and the observation will be greater, 
though the core radii of most galaxies are thought to be small enough to make SIS a good 
approximation. 



S 



The greatest difficulty in this or similar work is the small number of observed lens systems 
suitable for this kind of statistical analysis. We need to have more "clean" lens cases, 
meaning a single galaxy acting as lens. We hope that projects like Sloan Digital Sky Survey 
will answer most of the problems we have now. 
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FIGURES 



FIG. 1. Probability distribution functions of image separations in multiple- image lens systems 
for various angular resolution biases: NO (solid line), HST (dotted line), PSF1.0 (short dashed 
line), PSF0.7 (long dashed line), EYE1.0 (short dash-dotted line), EYE0.7 (long dash-dotted line), 
and MG-C (long dash-short dashed line). Larger triangles mark the image separations in the lens 
systems PG1115, Q2237, Q0142, H1413, MG0414, and B1422. Smaller triangles mark those in 
Q1208 and MG2016. 

FIG. 2. Cumulative distribution functions (CDF's) — fraction of cases with separations smaller 
than A6 — for various resolution biases and those from observed lens cases: (a) the dot-dashed curve 
is CDF for no selection bias (NO), the solid line CDF from lens systems O-ALL, the dotted line from 
O-SINGLE, and the dashed line from O/R-SINGLE; (b) the dotted curve for bias HST and the solid 
line from corresponding lens systems; (c) the dotted curve for bias PSF1.0, the short-dashed curve 
for bias PSF0.7, the long-dashed curve for bias EYE1.0, the dot-dashed curve for bias EYE0.7, and 
the solid line from corresponding lens systems; (d) the dotted curve for bias MG-C and the solid 
line from corresponding lens systems. 

FIG. 3. PDF's with [solid] and without [dotted] the (3/2) 1 / 2 correction factor for no selection 
bias (a) and for MG-C selection bias (b). 
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TABLES 



TABLE I. Relevant multiple image gravitational lens systems. 



Lens system 


A6» a 


zs 


PG1115 b 


2.3 


1.72 


Q2237 b 


1.8 


1.69 


Q0142 b 


2.2 


2.72 


H1413 b 


1.1 


2.55 


MG2016 C 


3.8 


3.27 


Q0957 d 


6.1 


1.41 


MG0414 e 


3.0 


2.64 


B1422 e 


1.3 


3.62 


Q1208 f 


.47 


3.80 


MG1131 S 


2.1 


1.13? 


MG1654S 


2.1 


1.74 



a Maximum image separation. 

b Discovered in an optical observation. Lensed by a single galaxy. 
c Discovered in an optical observation. Lensed by two galaxies. 
d Discovered in an optical observation. Lensed by a galaxy plus cluster. 
c Discovered in a radio observation. Lensed by a galaxy. 
Classified ELS cl proposed case fig] . 
g Radio ring. 
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TABLE II. Angular selection bias and corresponding lens systems. 



Selection Bias 



Lens Systems 



NO a 


1115, 


2237, 


0142, 


1413, 2016, 0957 (0-ALL b ) 


NO a 


1115, 


2237, 


0142, 


1413 (0-SINGLE c ) 


NO a 


1115, 


2237, 


0142, 


1413, 0414, 1422 (0/R-SINGLE d ) 


HST C 


1115, 


2237, 


0142, 


1413, 1208 


PSF1.0 f 


1115, 


2237, 


0142, 


1413 


PSF0.7 g 


1115, 


2237, 


0142, 


1413 


EYE1.0 h 


1115, 


2237, 


0142, 


1413 


PSFO.? 


1115, 


2237, 


0142, 


1413 


MG-C j 


2016, 


0414, 


1131, 


1654 



a No selection bias. 

b Lens systems discovered in optical observations. 

c Lens systems discovered in optical observations. Lensed by a single galaxy or lens not known. 
d Lens systems discovered in optical or radio observations. Lensed by a single galaxy or lens not 
known. 

c Selection bias for the HST snapshot survey. 

Selection bias for ground-based optical observations with the PSF subtraction method and a seeing 



FWHM=1.0" p|. 

g Selection bias for ground-based optical observations with the PSF subtraction method and a seeing 
FWHM=0.7" p|. 



h Selection bias for ground-based optical observations with visual examination and a seeing 
FWHM=1.0" p|. 



Selection bias for ground-based optical observations with visual examination and a seeing 

FWHM=0.7" g2|. 

j MG C band survey ll|J| 



TABLE III. Monte-Carlo simulation results. 



Selection Bias/Samples 


With (3/2) 1 /2 


Probability a 


Without (3/2) Va 


NO-O-ALL 


.01 




.00 


NO-O-SINGLE 


.37 




.04 


NO-O/R-SINGLE 


.30 




.01 


HST 


.56 




.11 


PSF1.0 


.47 




.08 


PSF0.7 


.42 




.06 


EYE1.0 


.52 




.12 


EYE0.7 


.44 




.07 


MG-C 


.09 




.00 


a Fraction of MC draws having likelihood smaller than that of the observed dat 


a. 
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TABLE IV. Kolmogorov-Smirnov test results. 



Selection Bias/Samples Confidence of Rejection 3, 

With (3/2)V2 Without (3/2) x l 2 



MO C\ ATT 


Q/l 


.yyo 


NO-O-SINGLE 


.71 


.97 


NO-O/R-SINGLE 


.87 


.995 


HST 


.17 


.81 


PSF1.0 


.53 


.92 


PSF0.7 


.60 


.94 


EYE1.0 


.39 


.89 


EYE0.7 


.50 


.92 


MG-C 


.98 


.999 



a Confidence to reject the null hypothesis that data come from the model distribution. 
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(a) 



TTT 



TT 



A 



I .J. 
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(C) 



- /// 



-r 
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(b) 



TTT 



TT 



- (d) 



I I 'I I I I I I I 



I I I I I I I I I I I I 



1 2340 1 234 

A0 (arcsec) 



